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Knudsen gas provides nanobubble stability 
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We provide a model for the remarkable stability of surface nanobubbles to bulk dissolution. The 
key to the solution is that the gas in a nanobubble is of Knudsen type. This leads to the generation 
of a bulk liquid flow which effectively forces the diffusive gas to remain local. Our model predicts 
the presence of a vertical water jet immediately above a nanobubble, with an estimated speed 
of ~ 3.3 m/s, in good agreement with our experimental atomic force microscopy measurement of 
~ 2.7 m/s. In addition, our model also predicts an upper bound for the size of nanobubbles, which 
is consistent with the available experimental data. 

PACS numbers: 



Classical diffusion predicts a lifetime of ^ 1 /is for a 
nanoscopically-sized bubble. So the fact that surface 
nanobubbles [l|-[lot (typical height ~ 20 nm and width 
^ 100 nm) persist for at least 11 orders of magnitude 
longer than this [llj is both remarkable and puzzling. 
Are classical diffusion laws simply not applicable at these 
length scales? Are nanobubbles coated with diffusion- 
limiting molecules [H, 13|? Or does the gas indeed dif- 
fuse out, but is balanced by an equivalent influx 
Supersaturation was thought to be the key to nanobub- 
ble nucleation and stability [isl . [IgI. but it is now known 
that this is not a requirement |17| . It was also origi- 
nally thought that nanobubbles may not actually contain 
gas [3] J but this is not correct either [l^ . This is one of 
the outstanding questions in fluid dynamics [13, H^l ■ 

Surface nanobubbles are fundamentally interesting. 
For example the gas-side contact angle that they make 
with the solid is always very low, regardless of the sub- 
strate chemistry, and also depends on size [2l| - [2^ . This 
is in clear contrast to the classical view that the contact 
angle is a material property and should be substrate de- 
pendent and size independent. Nanobubbles also have 
clear potential in applications, such as controlling slip 
in microfluidie devices l25l. \2& an d surface cleaning in 
nanofabrication processes [271, l28| . Hence, understand- 
ing their stability is paramount. 

In this Letter we suggest a solution to the mystery 
of nanobubble stability by demonstrating that the ex- 
act nature of the gas, i.e. Knudsen, is the key. The 
symmetry-broken geometry provided by the hard sub- 
strate and the 'leaky' liquid/gas interface thus generates 
a bulk gas flow. In turn, due to the continuity of shear 
stress boundary condition, this bulk gas flow leads to a 
bulk liquid flow and, due to conservation of mass, the 
gas-rich liquid circulates from the bubble apex around to 
the three-phase line. Hence, although the gas molecules 
do indeed dissolve into the liquid as expected, they re- 
main local to the bubble in the gas-rich liquid stream 
and are effectively transported back to the three-phase 
line for re-entry. We validate our theory by performing 



non-contact-mode open-loop atomic force microscopy in 
the liquid environment, which has enabled us to measure 
an incredible 2.7 m/s upthrusting water jet immediately 
above a nanobubble, in good agreement with the 3.3 m/s 
jet predicted by our model. Hence, our measurements 
clearly demonstrate that surface nanobubbles are in a 
dynamic equilibrium. 

Knudsen gases differ from their classical counterparts 
insomuch as the molecules hardly interact with each 
other. Thus, rather than a test volume possessing zero 
preferred direction, as would be the case for an ideal gas, 
Knudsen gas molecules mainly travel due to energy ex- 
change with the walls, so are heavily dependent on the 
geometry of their surroundings. Thus, the symmetry- 
broken geometry offered by the hard substrate and liq- 
uid/gas interface of a Knudsen gas-filled surface bubble 
leads to more gas traveling away from the solid substrate 
than is reflected back from the liquid/gas interface (i.e. 
some gas molecules are transmitted and diffuse away). 
Hence, the gas flow mimics the broken symmetry of the 
geometry such that every volume element of gas possesses 
a bulk flow away from the substrate. This is a generic 
statement for any Knudsen gas with thermal drive and 
one leaky wall (for nanobubbles the thermal drive comes 
from the substrate, which we treat here as a heat bath). 

The requirement for Knudsen gas behaviour is depen- 
dent on the Knudsen number, Kn, i.e. the ratio between 
the molecular mean free path. A, and the typical length 
scale of the container. For a surface bubble, this length 
scale is the bubble height, h, so Knudsen gas behaviour 
is exhibited when 
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where kT is the thermal energy, R the radius of curva- 
ture, 7 the surface tension, cr the molecular coUisional 
cross-section, pq the ambient pressure, and 9 the gas-side 
contact angle. For macroscopic bubbles R ^ 2'y/po so 
Kn ^ 1 and the gas is ideal, but for nanoscopic bub- 
bles R <^ 27/po and the requirement for Knudsen gas 
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behaviour becomes surprisingly only dependent on the 
contact angle. As an example, an oxygen nanobubble on 
graphite at room temperature will exhibit Knudsen gas 
behaviour if the contact angle is ^ 25°. 

Coincidentally, an as yet unexplained mystery of sur- 
face nanobubbles is that their contact angles are not con- 
sistent with their micro- or macroscopic counterparts. In- 
stead, their contact angles are always found in the range 
5°^6'^25° for hy drop hobic surfaces, regardless of the 
substrate chemistry |2ll424l |. suggesting that their inter- 
nal gas is always of Knudsen type. In what follows we 
shall demonstrate that this puzzling observation may be 
the key to the nanobubbles' mysterious stability. 

We begin with the gas molecules arriving at the liq- 
uid/gas interface. These molecules have not been able to 
interact with each other, due to their Knudsen behaviour, 
so they still possess the directional flow away from the 
substrate (Figure [T^). At the apex of the nanobubble 
they either reflect back into the bubble or transmit and 
dissolve. However, everywhere else on the liquid/gas in- 
terface the gas molecules arrive, on average, non-normal, 
so their bulk velocity can be decomposed into both a 
normal and a tangential component. The normal com- 
ponent acts like at the apex (reflection or transmission), 
but it is the tangential component that we are interested 
in. If a tangential velocity component exists in the gas 
phase it always points towards the apex of the bubble. 
This then communicates with the liquid phase through 
the assumed continuity of shear-stress boundary condi- 
tion (Figure [TJd), i.e. 



^J-g 



du 



du 
dn 





(2) 



FIG. 1: Knudsen gas streaming and nanobubble geometry, 
(a) The broken symmetry created between one solid surface 
and one 'leaky' liquid/gas interface leads to bulk upwards flow 
in the Knudsen gas. (b) The tangential component of the bulk 
Knudsen gas flow drives a bulk liquid flow due to the conti- 
nuity of shear stress at the liquid/gas interface, (c) Finally, 
due to conservation of mass, this gas-rich liquid stream cir- 
culates upwards at the bubble apex and back around to the 
three-phase line, effectively transporting the diffusive outflux- 
ing gas back to the three-phase line for re-entry p^ . 



where is the viscosity and ^ is the tangential velocity 
gradient of the gas, 5, or liquid, phase. In this pic- 
ture we have treated the system as a continuum (we are 
using Kn = 1, which means the gas is 50% ideal and 
50% molecular, both interpretations are correct) - this 
macroscopic assumption is at least partially valid but it 
must eventually breakdown. (The microscopic descrip- 
tion is complimentary: The gas molecules in the bubble 
have a net upflux, but they must travel radially outwards 
once they dissolve in the liquid. In order to change direc- 
tion at the liquid/gas interface they transfer momentum 
to the liquid molecules in the direction from the three- 
phase line towards the apex.) Hence, the upward flow 
of the Knudsen gas in the nanobubble induces a bulk 
flow in the adjacent liquid. The strength of the drive is 
dependent on the precise position on the liquid/gas inter- 
face, with maximum effect near the three-phase line and 
zero effect at the apex. The importance of this result is 
shown in Figure [It: The flow in the liquid tangential to 
the liquid/gas interface creates a local circulatory stream 
due to mass conservation. (Note that diffusive gas flow 
can also lead to liquid streaming in the extreme limit of 



inviscid flow [29|.) 

The appropriate scaling of the shear stress is the molec- 
ular speed distributed over the bubble length scale, i.e. 

~ Ug/R, whilst in the liquid the induced velocity 
is distributed over the radius, L, of the circulation loop, 
i.e. f^l; ui/L. So, although the gas in a nanobub- 
ble does indeed dissolve into the liquid, this gas-rich liq- 
uid is streamed at speed u/ ~ from the bubble 
apex around to the three-phase line, where it can re-enter 
through either the attractive potential of the hydropho- 
bic wall [131 or through adsorption to the substrate and 
surface diffusion. 

To test our hypothesis of gas-driven streaming in the 
bulk liquid we have performed non-contact-mode open- 
loop feedback-disabled atomic force microscopy in the 
liquid environment. This has allowed us to make a di- 
rect measurement of the force-field in the liquid, with the 
AFM cantilever acting as a local force, and thus velocity, 
probe. The liquid was ultrapure water (Simplicity 185 
system, Millipore, France), which had been thoroughly 
degassed before supersaturation with 3 atm of argon gas. 



3 




This was then deposited onto freshly cleaved HOPG be- 
fore scanning in both tapping mode and the non-contact 
open-loop feedback-disabled mode. The cantilevers were 
Au-back-coated Si3N4 Veeco NPG probes (radius of cur- 
vature 30 nm, full tip cone angle 35°), with resonance 
frequencies in liquid of oJq' w 15 — 25 kHz. 

The results of our measurements are shown in Figure 
where the blue curve is a line scan over the apex of the 
nanobubble in Figure[2ji, and the red curves are measures 
of the force field taken at various distances above the 
substrate (from ^ 250 nm to ^ 600 nm, in steps of 
50 nm). A clear upthrust exists in the immediate vicinity 
of the nanobubble. The maximum force exerted by this 
fiow was found to be ~ 1 nN at a distance of 250 nm 
above the nanobubble, with the flow still measurable as 
far away as ^ 500 nm. 

(To test whether the upthrust force was electrostatic, 
we changed the bias of the cantilever with respect to the 
sample from —1.5 V to +1.5 V in steps of 0.5 V. However, 
no change in force was evident. We also investigated the 
effect of different ionic concentrations within the liquid by 
using 0.00 M, 0.01 M, 0.02 M, 0.05 M, 0.10 M, 0.20 M, and 
0.50 M NaCl solutions. Again, no change in the upthrust 
force was detectable. Thus, the upthrust forces measured 
from nanobubbles are not electrostatic.) 

The implication is clear: Our experimental observa- 
tions demonstrate that nanobubbles are not in a true 
equilibrium. This settles the ongoing debate in the 
field 



12l |14| : The system is in a dynamic equilibrium. 



We benefitted greatly in our experiments from the fact 
that the upthrusting jet was focussed immediately above 
the nanobubble. Due to the continuity of mass, a central 
up-flow must be balanced with a downwards flow in a 
circulatory stream. However, we were not able to mea- 
sure this downwards motion of the circulation stream, 
presumably because it is averaged out over a very large 
annular area. (The prediction from the continuity of flux 
is that the downwards velocity is ~ 1 % that of the up- 
wards jet, where wc have taken a fixed radius circulatory 
stream of 50 nm rising upwards immediately above the 
nanobubble and returning back downwards at a radial 
position of 2L.) 

We can validate our model further by measuring the 
velocity of the jetting water. Treating the nose of the 
AFM cantilever as a sphere of radius r, the effective fiuid 
velocity can be estimated by equating our 1 nN mea- 
surement with Stokes' drag, i.e. ui = F/Qnfir. In the 
case of the nanobubble in Figure [21 the water jet stream 
induced by the nanobubble's Knudsen gas travelled at 
an incredible 2.7 m/s! This is in good agreement with 
the ui ^ = 3.3 m/s prediction of our model, where 

we have used L = 600 nm — 90 nm as the radius of the 
circulation stream (the limit of the circulation loop mea- 
sured in Figure [5] was 600 nm above the substrate, and 
we have subtracted the 90 nm height of the nanobubble 
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FIG. 2: (a) Three-dimensional topological image (1 /im x 
1 X 100 nm) of one of the nanobubbles investigated in 
this study. The nanobubble was created by depositing room 
temperature ultra pure (Millipore) water on highly-oriented- 
pyrolitic-graphite (HOPG), also at room temperature. The 
water had been thoroughly degassed before supersaturating 
it with 300 % argon, resulting in a bubble with radius of 
curvature R ~ 1.4 /^m and contact angle 6 ~ 20°. The 
nanobubble was topologically imaged regularly over ~ 12 hrs 
to confirm that its size was not changing with time, (b) 
Cross-sectional line scan (blue) of the nanobubble in (a), with 
feedback-disabled non-contact-mode force field measurements 
(red), taken at 250 nm to 600 nm above the substrate, in 
steps of 50 nm. At 250 nm from the substrate (~ 160 nm 
above the nanobubble's apex) the AFM cantilever was de- 
flected upwards by a '-^ 1 nN force, which decreased in mag- 
nitude with increasing cantilever-substrate separation until it 
became smaller than the AFM resolution (~ 5 pN) at 550 nm. 
The force ordinate on the right-hand side of the image is 
for the bottom-most force field measurement. The proceed- 
ing force field measurements have the same scale, but with 
their zeroes offset according to the height at which they were 
recorded (left-hand ordinate). We recalibrated the AFM can- 
tilever to the substrate between each force field measurement 
in order to maintain accurate separation with minimal drift 
for each scan. 
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- the precise limit is expected to be between L and 2L). 

We must clearly consider heat generation from this ex- 
ceptional jet. Heat generation is through viscous dissi- 
pation within the liquid, and has rate €dis 
Hence, the rate of heating is 



dn\l- 



AT/ At = — 7- 
cp on 



(3) 



which is an incredible lO^K/s (where p and c are the 
water density and specific heat capacity, respectively). 
However, we clearly do not see such a large increase in 
temperature in the experiments. There are two possi- 
ble explanations: (i) Firstly, the heat is generated in the 
small volume of liquid that forms the recirculation loop, 
but this is advected away through the entire droplet. 
This reduces the temperature increase by the ratio of 
these two volumes, i.e. approximately 10"^'^ K in our ex- 
periment, (ii) Secondly, the relevant time scale for the 
rate of heating is the travel time of one loop of the cir- 
culatory stream, t = 2ttL/ui. Given that we are treating 
our solid surface as a heat bath, and that the solid surface 
is much more efficient at conducting the heat away than 
the liquid, once the heated liquid has circulated around 
to the wall it can efficiently exchange energy and return 
to temperature T. The estimated temperature increase 
of the liquid using this time scale is sa 10 mK. In either 
case, dissipative heating is negligible. 

Finally, our model predicts an upper limit for the size 
of a nanobubble. We have chosen to use Kn = 1 as the 
limiting factor for the stabilising effect in this Letter. If 
the height of the nanobubble was larger than the mean 
free path of gas in atmospheric conditions (Ag ~ 100 nm), 
we may expect the Knudscn gas behaviour to break 
down. All nanobubble studies to date have had nanobub- 
bles with heights smaller than this. 

To summarise: (i) Surface nanobubbles contain Knud- 
sen gas which possesses a bulk flow due to the thermal 
energy of the substrate and the leaky liquid/gas interface; 
(ii) this bulk gas flow drives a bulk liquid flow due to the 
continuity of shear stress; (iii) the gas-rich liquid flow is a 
circulatory stream from the apex to the three-phase line, 
due to the conservation of mass; (iv) the gas re-enters 
the nanobubble at the three-phase line, replenishing the 
diffusive outflux. 
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